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It is shown that a nonclassicality witness, whose expectation value can be measured for all quantum states, can 
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I. INTRODUCTION 

The notion of nonclassicality describes the difference be- 
tween classical and quantum physics. For the harmonic oscil- 
lator, its formal definition is based on the coherent states \a), 
which are the closest analogues to the classical oscillation JT]. 
Any quantum state p can be formally written as a statistical 
mixture of (classical) coherent states 0O, 

p = J d 2 aP(a)\a)(a\. (1) 

The function P(a) is the so-called Glauber Sudarshan P func- 
tion. If it is a classical probability density, the state p is a 
classical mixture of coherent states. In many cases, however, 
the P function is a quasiprobability, violating properties of a 
probability. Then the state is called nonclassical J4|. 

In general this definition cannot be applied in most practical 
situations, since the P function may have strong singularities. 
Therefore, several criteria for nonclassicality have been de- 
veloped, which can be distinguished in two major categories. 
First, one can examine the expectation value of a so-called 
nonclassicality witness W for a given quantum state |3][6]. 
Without loss of generality, this operator has a nonnegative ex- 
pectation value for all classical states p c \, 

Tr{ Pcl W} > 0. (2) 

Conversely, if this inequality is violated for a quantum state p, 
then this state must be a nonclassical one: 

Tr{pW} < => p nonclassical. (3) 

Based on nonclassicality witnesses, one may formulate condi- 
tions for matrices of characteristic functions [7|, moments |8] 
9 1, or outcome probabilities iflOl . 

On the other hand, one may examine phase-space repre- 
sentations of a quantum state, which play the role of a prob- 
ability density in some sense. Besides the P function, the 
Wigner function is frequently used in experiments since it is 
always regular. For nonclassical states, the Wigner function 
may become negative for some points in phase space. These 
negativities indicate the nonclassicality of the quantum state. 
However, nonnegativity of the Wigner function does not prove 
classicality, and negative values of the Wigner function are 
only sufficient for nonclassicality. The recently developed 



concept of nonclassicality filters and quasiprobabilities pro- 
vides a complete characterization of nonclassical effects ifTTTl 
in terms of regular phase-space distributions. 

In the present paper we show that generally applicable non- 
classicality witnesses can be constructed from any nonclas- 
sicality filter. In Sec. II, we elaborate the connection to the 
previously known nonclassicality witnesses and provide an 
explicit analytical expression of a universal witness. Exper- 
imental applications are briefly discussed in Sec. III. In Sec. 
IV, the usefulness of our method is illustrated for an example 
of a quantum state, under conditions when many other tests 
fail. 



II. UNIVERSAL NONCLASSICALITY WITNESSES 

A. Nonclassicality filters and witnesses 

Let us start from the characteristic function of the P func- 
tion of the state p, 

*(/?) = Tr{pe /5a V ra } . (4) 

We choose a function of a complex argument (5, £l„ (/3) = 
fii(/3/w), satisfying the conditions ifTTI : 

CI. ri w (/?)e' /3 ' / 2 is integrable for all positive w. 

C2. O w (/3) has a nonnegative Fourier transform. 

C3. Q w (0) = 1 and lim,,^ n w {(3) = 1, for all /3. 

Then, f2 w (/3) is a family of nonclassicality filters, parameter- 
ized by the real number w. 

The Fourier transform of the filtered characteristic function, 

Pw(a) = ^ J d 2 /3$(/3)a,(/3)e^*^ Q * /3 . (5) 

is a filtered P function of the state. Since the filter is designed 
for the requirements of nonclassicality tests, P w (a) is denoted 
as nonclassicality-filtered P function (NFP). The NFP is non- 
negative for all classical states. For all nonclassical states, 
its negativities uncover nonclassicality for properly chosen w 
and a. When the NFPs contain the full information about the 
quantum state, they are called nonclassicality quasiprobabili- 
ties IfTTI . Their experimental reconstruction has been demon- 
strated lfT2l[l3l . In the following we renounce this additional 
constraint and consider general NFPs. 
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B. Construction of nonclassicality witnesses 

In the following, we rewrite the NFP as the expectation 
value of an Hermitian operator with finite trace. For this pur- 
pose, we insert Eq. Q into the definition of the NFP Q and 
obtain 

P w (a) = Tr [pW w (a)} = (W w (a)), (6) 
with the operator 

W w (a) = ^ J d 2 /3ft l ,(/3)e /3(at - Q * ) e- /r ( a - Q >. (7) 

The complex argument a plays the role of a coherent displace- 
ment of the operator W w = W w (0), it can be expressed with 
the unitary displacement operator D(a) = e aa ~ a a by 

W w (a) = D(a)W w D(-a), (8) 

Let us consider the operator W w first. Because £1 W (J3) has a 
real Fourier transform, we have £l vv (— (3) = fi*(/3), which 
implies that W w is Hermitian. Furthermore, the expectation 
value for any coherent state is given by the Fourier transform 

(a\W w \a) = ^ J d?p n w (py a *-P" a , (9) 

which is nonnegative by definition of a nonclassicality filter. 
In consequence, the same holds not only for coherent states, 
but for all classical states. Last, W w (a) is a trace class opera- 
tor, since Tr{e /3at e- /3 * a } = tt<5(/3) and fi w (0) = 1: 

Tr{W w } = — 2 /d 2 /3a,(^)Tr{e^ t e- /3 * a } = -. (10) 

Since the coherent displacement does not affect these proper- 
ties, we conclude that W w (a) has the same features for all a. 

This gives rise to the interpretation of the Hermitian opera- 
tor W w (a) as a nonclassicality witness: Its expectation value 
is nonnegative for all classical states. Hence, states with nega- 
tive expectation values are clearly identified as being nonclas- 
sical. In this sense, we have proved that any nonclassicality fil- 
ter £l w (0) can be used to define a nonclassical witness W^fa) 
by Eq. n\. Its expectation value is the NFP P w (a), see Eq. d6). 
Finally, by adapting the results in ifTTl . we know that for any 
nonclassical state and any nonclassicality witness W w (a) as 
defined above, there exists a real width parameter w and a 
complex number a such that the expectation value of W w (a) 
is negative. Therefore, we constructed a set of nonclassicality 
witnesses, parameterized by only three real numbers, which 
is sufficient for the detection of nonclassicality of an arbitrary 
quantum state. In this sense, any witness constructed in this 
way is universal. 

In order to demonstrate the equivalence of witnesses and 
filters, we still have to show the converse statement: When- 
ever one has a nonclassicality witness W with Tr(W^) = —, 



as will be assumed in the following, one may construct a non- 
classicality filter from it. For this purpose, let us represent 
the witness W by its characteristic function <fr^ (/?) of the Q 
function JTU, 

W=~ [ d 2 l3<5> { 9\/3)e^e- ' & . (11) 

TT Z J W 

Now, one easily sees from Eq. |7]) that the characteristic func- 
tion can be directly connected to a nonclassicality filter: 

*U8) = *S?G8)> d2) 

and the corresponding NFP is given by Eq. Q. Because of 
its bounded trace, the expectation value of W is always finite. 
Hence the NFP is always a regular function, and may be ac- 
cessible to experiments. Still, the width parameter w is not 
implemented, since the right hand side of Eq. (JT2J does not 
depend on a real parameter w. However, in Appendix A we 
show how this can be achieved. 



C. Explicit form of universal witness operators 

Let us consider the nonnegativity of the Fourier transform 
of the filter Vl w (f3) more in detail. This condition can by writ- 
ten with the help of Eq. Q as 

(a\W w \a) =til(a*,a)& w (a*,a), (13) 

with some suitable function Q w (a* ) a). Here, we write the 
arguments a and a* explicitly in order to have a clear notation 
of corresponding operators which will be defined below. The 
nonclassicality filter is the autocorrelation function 

n K (/3) = J d 2 p'u w (p')uj w 08 + tf) (14) 

of the Fourier transform of to w (a* , a), 

"wifl = \J d 2 aQ w {a*,a)e a ^- a ^' . (15) 

Note that Q w (0) = 1 sets a normalization condition on 
w w (a* , a). However, it is not needed for nonclassicality tests, 
but only for the correct normalization of quasiprobabilities. 
Furthermore, since we consider fi w (/3) = Qi(/3/w), the width 
parameter enters into ui(/3) as 

uj w (a* , a) — wuji(wol*, wa). (16) 

Most importantly, u){a*,a) has to be chosen such that the 
function f2 w (/J) still fulfills condition CI. It has been shown 
to be sufficient that u) w (/3) itself fulfils this requirement IfTTl . 
i.e. w„,(/3)e'' 3 ' I" 1 is integrable for all w. 

From Eq. (13) , we easily find the normally-ordered form of 
the nonclassicality witness W w . It is obtained by replacing the 
c-numbers a, a* by the operators a, a): 

W w =:Ql(at,a)u> w (al,a) :, (17) 
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the notation : : indicates that all products of the operators a) , a 
have to be written in normal order. Furthermore, by defining 
the displaced operator 



u>w,a(or, a) — u) w (w — a* ,a — a), 
the relation ([8} leads us to 

W w (a) =:uil a (a t , a)w M ,, a (d t , a) : 



(18) 



(19) 



This is a remarkable result. First, we can construct nonclas- 
sicality witnesses easily by properly choosing the function 
Q w (a*,a). It has to satisfy some weak conditions, in par- 
ticular its Fourier transform obeys the condition CI. The 
other two conditions are then trivially satisfied by construc- 
tion. Second, it states that any nonclassicality witness, which 
can be constructed by a nonclassicality filter, can be writ- 
ten in such a form. The operator Q wa (a',a) can be found 
from Eqs. (|9j, ( fT3| , ( fT8| ), and only depends on three real 
parameters. As already discussed, one particular operator 
£y a (dT,d), equipped with these three parameters, is suffi- 
cient to completely determine nonclassicality of an arbitrary 
state. 



D. Example of a universal witness 

Here we consider a witness that can be given explicitly. We 
fix w vv (a^,d) by choosing its Fourier transform as the disc 
function 



1, |j3| < 1/2, 
0, elsewhere. 



(20) 



The corresponding nonclassicality filter does not preserve all 
information about a quantum state, since it has bounded sup- 
port. Furthermore, it is not appropriately normalized. How- 
ever, it is suitable for the construction of a universal set of 
nonclassicality witnesses. 

From the inverse relation of ( flB) and Eqs. ([16), ( p"8j ), (jT9j, 
we find the nonclassicality witness: 



W w (a) =: 



WWW -a*)(o-q))] s 
4(dt — a*)(d — a) 



(21) 



where J\ (x) is the Bessel function of first order. Equation ( |2"T| > 
represents a complete family of witness operators in a closed 
analytic form, which is an important finding of our paper. It is 
suited to identify all the nonclassical effects of any quantum 
state of a harmonic oscillator. This expression is well-defined 
in general. It can be expanded into a normally ordered power 
series of the displaced photon number operator n(a) = (a* — 
a*)(a — a), 



W w {a) 



w 



(-w 2 /4) m (2m + 2 



16 ^ [(m + l)!l 2 

m— 



: h( a y 



(22) 



It is noteworthy that this form of the witness is not unique. As 
we have shown, different witnesses can be constructed from 
different nonclassicality filters, defined by different functions 
Wi(/3). All these witnesses are equivalent for the task to un- 
cover nonclassical effects completely. 



E. Relation to other witnesses 

Nonclassicality witnesses of a form similar to (jT9]» have al- 
ready been examined in [5,6, 8|. They have been constructed 
as normally ordered squares of some operator / = f(cr,a), 

w =:/t(at,a)/(at,a):. (23) 

Using the P function, its expectation value is given by 

(:p(a\a)f(a\a):)= f P(a)\f(a*,a)\d 2 a. (24) 



Clearly, negative expectation values can only be achieved by 
negativities of the P function, i.e. by nonclassical states. Sim- 
ple examples lead to well known nonclassicality conditions. 
Choosing / = x — (x), where x is the quadrature operator, we 
may identify squeezing 1151 . For / = h — (n), h = a* a, we 
uncover sub-Poissonian photon number statistics ifTBl . 

More general nonclassicality tests can be obtained from 



Eq. (23 1. One may derive a complete hierarchy of conditions 
in terms of characteristic functions [7 1, which directly follows 
via Fourier representation of the operator /(d^,d), cf. [8|. 
Via Taylor series expansion of /(d^,d) one obtains criteria 
in terms of matrices of moments [5 , 8 9 1. The former method 
can be applied to any quantum state. In general, however, such 
a nonclassicality test may become rather complex. The latter 
method is only applicable provided that the needed moments 
exist. It also becomes complex, in particular if large matrices 
are needed to detect the nonclassical effects. If all moments of 
a quantum state exist, the test operators / can be restricted to 
polynomials in a, d^ [6 1, which, however, still form an infinite 
dimensional space. 

In contrast, our approach is much simpler, since it requires 
only control of three real parameters, in order to perform 
complete nonclassicality tests. They enter into the opera- 
tor uv Q (a^, d) simply as a displacement and rescaling. On 
the other hand, our method does not rely on the assumption 
that all moments of the quantum state exist. This is due to 
the requirement CI for the nonclassicality filter, which guar- 
antees the existence of the quasiprobabilities and expecta- 
tion values of the witnesses. Remarkably, this condition is 
much more fundamental as it might look at first sight. It is 
also necessary for a complete test of nonclassicality of states, 
for which all moments exist. For instance, if we choose 
/w,a(d^,d) = w(d — a), it has the same structure as our 
w, v ,q (d^ , d) in Eq. dTSJ, but cannot be used to construct a non- 
classicality filter which fulfils the condition CI. As a witness, 
with / = f w<a in Eq. (23 1, it requires the existence of the first 
and second moments of d^, a. However, this does not provide 
a nonclassicality test, even for a single quantum state! 



III. DIRECT MEASUREMENT OF NONCLASSICALITY 
WITNESSES 

From Eq. ([6|, it is clear that one may calculate the expec- 
tation value of the witness W w (a) by reconstruction of the 
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NFP P w (a), analogously to the work in lfl2l [131 . However, 
this procedure requires full quantum tomography of the state. 
The question arises if one can estimate the expectation value 
directly. 

Let us restrict to the case of phase-independent nonclassi- 
cality filters, i.e. O w (/3) = f2 u .(|/3|). The resulting nonclassi- 
cality witness W w is diagonal in photon number basis, cf. e.g. 



Displacement 



Detection 



Eq. (22 1. Its expectation value can be written as 



oo 

E 

n=0 



(W w ) = J2(n\W w \n)p n , 



(25) 



with p n being the photon number statistics of the quantum 
state. Since the latter can be measured by photon-number re- 
solving detectors [17], and the matrix elements of W w can be 
calculated theoretically for arbitrary w, one can determine this 
expectation value with standard experimental techniques. 

Of course, practically one is restricted to a finite number 
of photons, approximating the series in Eq. ( |2"5j ) by a finite 
sum. The smaller the width parameter, the better the approxi- 
mation will be, since the matrix elements (n|Wlv|n) are poly- 
nomials in w 2 , cf. Eq. ( f22) . For larger w parameter, a larger 
number of photons must be resolved, which can be difficult in 
practice. This problem is similar to the practical limitation of 
the w parameter as discussed in the context of nonclassical- 
ity quasiprobabilities [12|. In such cases one may estimate a 
systematic error due to the truncation of the series d25). 

In order to measure the expectation value of W H .(a), we still 
have to include the coherent displacement in Eq. dS). For this 
purpose, it is indifferent if we displace the witness operator by 
the amplitude a or the quantum state by the amplitude —a, 



(W w (a)) 



OO 

£< 

n=0 



n\W w \n)p n (-a), 



(26) 



with p n (a) = (n\p a \n) and p a = D(—a)pD(a). This can 
be done by methods considered in lfT71 - [T9l : a highly trans- 
parent beamsplitter is used to displace the original state by 
a coherent one and a photon-number resolving detector mea- 
sures the signal, see Fig.[T] Then, the expectation value of the 
phase-independent nonclassicality witness W w is determined 
for the displaced state, which equals to the expectation value 
of W w (a) for the original state. In this way, one can choose 
a by the experimental setup, while the width parameter is in- 
cluded in the Fock matrix elements of W w . The variation of 
these three real parameters enables one to perform a complete 
nonclassicality test for an arbitrary quantum state. Note that 
the practical determination of (W w (a)) can also be based on 
direct sampling of the photon number statistics from phase- 
randomized quadrature measurements, as done in ifTTl . 



IV. APPLICATION 

Let us consider a single photon added thermal state 
(SPATS) 1201 . characterized by the mean photon number of 
the thermal background, n, and a quantum efficiency rj. It 
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Figure 1. Experimental setup. The beamsplitter displaces the initial 
state, followed by a photon-number resolving detection. The filter 
width is introduced by application of Eq. §26). 
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Figure 2. (Color online.) Dependence of the expectation value of 
the witness W„,(0) on the filter width, for SPATSs of different mean 
thermal photon number n and quantum efficiency r\ — 0.5. Nonclas- 
sicality is indicated by negative expectation values, which appear for 
sufficiently large w for all states. 



is known that for n > l/\/2, this state does exhibit neither 



squeezing nor sub-Poissonian statistics, arising from Eq. (23 1 
with f = x — (x) and f — h — (n) respectively [5|. Fur- 
thermore, for n > 0.39, the characteristic function does not 
show nonclassicality of first order [7|. Therefore, these well- 
known nonclassicality witnesses are not able to verify non- 
classicality of this particular state. It is also important that for 
r\ < 1/2, the Wigner function is positive semidefinite ||2~TI . 
Hence, by setting the quantum efficiency r\ = 1/2 and choos- 
ing n > 1/V2, it is a nontrivial task to witness the nonclassi- 
cality of the state under study. 

For our examinations, we choose the nonclassicality filter 
given in Eq. ( p2| . Since SPATS are phase symmetric and sim- 
ilar to the photon, we may test the witness at the origin of 
phase space, i.e. a = 0. Figure [2] shows the dependence of 
the expectation value of the nonclassicality witness W w on 
the filter width w for states with different mean thermal pho- 
ton number n. It is clearly seen that for sufficiently large w, 
nonclassicality of all states is verified by negativity of the ex- 
pectation value of the witness. The lower the mean thermal 
photon number, i.e. the closer the state is to the photon, the 
more pronounced is the effect. 
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V. CONCLUSIONS 

We provide universal nonclassicality witnesses, which 
identify the quantum effects of any quantum state of the har- 
monic oscillator. Our witness operators are controlled by only 
three real parameters, compared to the dependence of previ- 
ously known witnesses on an infinite number of parameters. 
The detection of a witness requires a controlled coherent dis- 
placement of the quantum state and the determination of the 
photon number statistics. We demonstrate the witnessing of a 
photon on a thermal background, which works very well even 
when other methods fail. 
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we need an additional nonclassicality filter £l' H ,(f3), which can 
be constructed by an autocorrelation function as proposed 
in ifTTI . Now, we look for suitable functions f(w), g(w) such 
that 

n w (/?) = f(w)*$ (f ) e -sHi^i 2 /2 + (i - f( w m> w (f3) 

(Al) 

is a nonclassicality filter and f2i(/3) = 

The function f(w) shall be a weight function, i.e. < 
/(w) < 1. Furthermore, for w = 1, we want to return to 
the original nonclassicality witness, cf. condition (b). This 
requires f(l) = 1 and g(l) = 0. Therefore, the function 
f(w) — exp(— (w — l) 2 ) is an appropriate choice. 

In order to have a nonnegative Fourier transform of Sl w (,3), 



it is sufficient that each summand in Eq. (Al i has this prop- 
erty. Since £l' w {0) is a nonclassicality filter itself, the state- 
ment is clear for the second term. Furthermore, the Fourier 
transform of is nonnegative, since 



Appendix A: Construction of NFPs from witnesses. 

In the main article, we have already shown that any non- 
classicality witness W with Tr(VF) = — can be used to define 
a nonclassicality filter fii(/3). Here, we construct a complete 
family of nonclassicality filters Sl w (/3), which can be used for 
a complete nonclassicality test. In order to do so, we have to 
find a filter, equipped with a parameter w, and with the fol- 
lowing requirements: 

(a) The properties C1-C3 for a nonclassicality filter are sat- 
isfied. 

(b) For w = 1, we want to return to the nonclassicality 



witness, 



We have already seen that (f3) fulfills two conditions for a 
nonclassicality filter: The corresponding NFP is always finite 
for all a, as we required by condition CI, and its negativities 
clearly indicate nonclassical effects, see condition C2. 

The task is to introduce a width parameter w to the function 
such that the constructed family of functions satis- 
fies all properties of a nonclassicality filter. For this purpose, 



$ ( £r ) (fi)e afl '- a ' fl <Pp = (a\W\a). 



(A2) 



equals to the expectation value of the nonclassicality witness 
for the classical coherent states. Therefore, it is sufficient that 
the remaining factor e sC ll 1 ^1 / 2 nas a nonnegative Fourier 
transform, which can be achieved by requiring g(w) > 0. 

Last, we require that fl w (/3)e^^ / 2 is integrable for any 
real w > 1, see condition CI. While the second term in 
Eq. ( |A1[ > fulfills this condition by its definition, the first term 
has to be examined more carefully. We already know that 
$|? ) (/3) e l' 3 l 2 / 2 is integrable, since fulfills condition 

CI. Therefore, we split the term as follows: 



W K w 



$(Q) (§\ e \P\ 2 /(2w 2 ) 



Jl-g(w)-l/w 2 )\f3\ 2 /2 



(A3) 



As the factor in square brackets is integrable, the remaining 
factor must be bounded by one, which means g(w) > 1 — 
1/w 2 . Altogether, the function g(w) = max{l — l/w 2 ,0} is 
a good choice, such that f2 u .(/3) defines a nonclassicality filter. 
This concludes our construction. 
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